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ABSTRACT 


A  resume  of  the  work  of  the  AV CO  Research  Laboratory 
during  essentially  the  first  year  of  its  existence  in  the  field  of 
high  temperature  gas  dynamics  is  presented.  This  laboratory  was 
organized  early  in  1955  by  Dr,  Arthur  Kantrowitz  to  investigate  the 
physical  gas  dynamic  problems  associated  with  the  high  tempera¬ 
tures  involved  in  flight  at  hypersonic  Mach  numbers  The  central 
concept  of  this  work  is  that  the  important  elements  of  the  environ¬ 
ment  anticipated  could  be  duplicated  in  shock  tubes 

The  advantages  of  the  straight  shock  tube  for  investigation 
of  high  temperature  gas  dynamics  are  discussed  For  many 
phenomena,  such  as  boundary  layer  problems  (other  than  transi¬ 
tion),  the  parameters  important  in  classical,  low  temperature 
aerodynamics,  such  as  Mach  number  and  Reynolds  number,  are 
shown  to  lose  much  of  their  importance  in  the  high  temperature, 
dissociating  gas,  situation  Instead,  the  flow  chemistry  is  shown 
to  be  a  critical  simulation  parameter  The  kinetics  of  the  flow 
chemistry  are  simulated  if  the  enthalpy  and  pressure  of  the  gas 
are  reproduced 

The  enthalpy  and  pressure  attainable  at  the  stagnation  point 
of  a  blunt  model  in  a  shock  tube  are  given  and  are  shown  to 
duplicate  those  encountered  at  the  stagnation  point  in  flight  up  to 
the  satellite  velocity  In  addition  the  condition  in  the  hot  gas 
behind  both  the  moving  shock  wave  and  behind  a  shock  wave  re¬ 
flected  from  the  closed  end  of  the  tube  are  shown  to  be  applicable 
to  the  study  of  many  hypersonic  flight  phenomena. 

Several  investigations,  conducted  in  shock  tubes,  into 
critical  gas  dynamic  problems  associated  with  the  high  tempera¬ 
ture  are  described.  Heat  transfer  measurements  in  the  laminar, 
dissociated,  boundary  layer  at  the  stagnation  point  are  discussed. 
Measurements  of  the  radiative  emissivity  and  electrical  conduc¬ 
tivity  of  high  temperature  air  are  described.  Other  investigations 
into  the  properties  of  high  temperature  gases,  including  prelimi¬ 
nary  studies  of  magnetohydrodynamics  are  summarized. 
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INTRODUCTION 


At  a  previous  meeting  of  the  AGARD  Wind  Tunnel  and  Model 
Testing  Panel  in  Paris  the  shock  tube  as  a  high  temperature  aero- 
dyn  amic  research  and  development  tool  was  discussed  by  Colonel 
Dodge^^^,  Much  of  the  discussion  centered  about  ways  of  producing 
hypersonic  Mach  numbers  to  completely  simulate  the  flight  environ¬ 
ment.  More  recently  in  Ottawa,  Dr  Kantrowitz^®^  addressed  this 
group  and  outlined  a  program  designed  to  answer  many  of  the  prob¬ 
lems  of  high  temperature  hypersonic  gas  dynamics  by  the  use  of 
straight  shock  tubes  without  high  Mach  number  nozzles,  and  within 
the  state  of  shock  tube  technology  as  it  then  existed 

Since  early  1955,  the  AVCO  Research  Laboratory  has  been 
working  on  implementing  this  program  under  the  direction  of  Dr, 
Kantrowitz  The  work  of  this  group,  over  essentially  the  first  year 
of  its  existence,  will  be  reported  in  this  paper  It  has  been  our 
philosophy  in  this  research  not  to  attempt  to  develop  a  facility  with 
the  maximum  flexibility  for  investigating  hypersonic  problems,  but 
to  consider  specific  areas  and  design  apparatus  which  could  be  quickly 
constructed  and  which  through  theoretical  analysis  and  experimental 
verification  may  shed  some  light  on  high  temperature  phenomena 

The  relative  kinetic  energy  of  the  air  intercepted  by  a  hyper¬ 
sonic  vehicle  flying  at  the  satellite  velocity  is  large  enough  to  dis¬ 
sociate  all  the  oxygen  molecules  into  atoms,  .dissociate  about  half  of 
the  nitrogen,  and  thermally  ionize  a  considerable  fraction  of  the  air 
The  presence  of  free  -atoms,  electrons ,  and  molecules  in  excited 
states  can  be  expected  to  complicate  heat  transfer  through  the  bound¬ 
ary  layer  by  additional  modes  of  energy  transport  such  as  atom  dif¬ 
fusion  carrying  the  energy  of  dissociation  Radiation  by  transition 
from  excited  energy  states  may  contribute  materially  to  radiative 
heat  transfer  There  is  also  a  possibility  of  heat  transfer  by  elec¬ 
trons  and  ions  The  existence  of  large  amounts  of  energy  in  any  of 
these  forms  will  undoubtedly  influence  the  familiar  flow  phenomena. 

Exploration  of  the  relationship  between  the  high  temperature 
phenomena  and  gas  dynamics  requires  careful  consideration  of  both 
the  environment  and  the  problem,  because  the  usual,  non  -  dis  sociating 
gas,  aerodynamic  parameters  are  not  necessarily  required  for  sim¬ 
ulation,  Aerodyttamic  heat  transfer,  like  other  boundary  layer 
problems,  does  not  require  the  simulation  of  Mach  number  ahead  of 
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the  bow  shock  so  long  as  the  pressure  distribution  is  similar  in 
the  test  and  flight  situations*.  This  can  easily  be  seen  by  con¬ 
sidering  the  boundary  conditions  for  the  boundary  layer.  In  the 
shock  tube  one  can  simulate  the  stagnation  enthalpy  and  pressure. 
Thus,  if  the  pressure  distribution  is  the  same  in  flight  and  shock 
tube,  the  absolute  values  of  the  gas  properties  at  every  point 
outside  the  boundary  layer  are  identical  for  the  two  cases  if  ther¬ 
modynamic  equilibrium  exists.  For  boundary  layer  problems 
simulation  is  thus  complete  in  equilibrium,  but  are  departures 
from  equilibrium  identical  ?  In  the  shock  tube  the  incident  shock 
may  leave  the  air  in  a  non-equilibrium  state,  i  e.  ,  partially  dis¬ 
sociated  However,  as  the  air  passes  through  the  bow  shock  at 
the  nose  of  a  model  it  is  immediately  raised  to  a  very  high  tem¬ 
perature  where  dissociation  proceeds  very  rapidly  The  high 
temperature  will  effectively  erase  any  memory  of  the  gas  state 
upstream  from  the  shock  wave,  since  dissociation  will  proceed 
very  rapidly  until  the  translational  temperature  is  reduced  to  a 
value  low  enough  to  slow  up  the  reaction  rate.  Thereafter  the 
chemical  state  of  the  gas  is  determined  by  the  pressure  time 
history  This  depends  only  on  the  initial  state  and  the  pressure 
distribution,  which  are  both  reproduced  in  the  shock  tube.  The 
chemical  state  is  therefore  simulated  correctly  at  every  point. 

However,  in  a  dissociated  gas  in  a  realistic  shock  tube, 
the  model  size  must  be  different  from  the  flight  vehicle,  and, 
because  the  degree  of  dissociation  depends  on  both  tempera¬ 
ture  and  pressure,  the  pressure  and  model  size  cannot  be  used 
to  compensate  each  other  as,  for  instance,  Reynolds  No.  simulation 
in  classical  aerodynamics  Thus  both  Reynolds  number  and  Mach 
number  lose  some  of  their  significance  as  s imulation  par ameter s . 

Realization  of  this  fact  eliminates  one  of  the  most  impor¬ 
tant  difficulties  in  the  application  of  the  shock  tube  to  hypersonic 
aerodynamics  Complete  simulation  is  theoretically  possible  in 
shock  tubes  by  the  use  of  the  expansion  nozzles  but  demandB 


*  It  will  be  shown  in  a  later  section  that  for  certain  blunt  bodies, 
such  as  a  hemisphere-  cylinder,  the  shock  tube  pressure  distribu¬ 
tion  is  Newtonian  and  therefore  identical  to  the  estimated  flight 
pressure  distribution  at  Mach  three  or  more  As  a  result*  th* 
pressure  distribution  simulation  is  independent  of  Mach  numbers. 
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excessive  driver  pressures  and  nozzle  expansion  ratios  of  the  order 
of  10^  However,  when  the  large  total  pressure  loss  across  a  high 
Mach  number  bow  shock  is  eliminated,  as  in  the  straight  shock  tube, 
the  smaller  degree  of  irreversibility  allows  stagnation  pressures  and 
pressure  distributions  typical  of  many  hypersonic  flight  situations  to 
be  achieved. 

For  radiation  studies  shock  tube  simulation  is  more  obvious 
It  has  been  stated  that  the  thermochemical  state  of  the  air  is  cor* 
rectly  reproduced  in  the  shock  tube.  Hov.;ever.  the  thickness  of  the 
hot  gas  region  is  much  smaller  and  the  conversion  to  flight  must  be 
made  by  geometric  consideration  of  the  two  cases 

Like  radiation,  other  physical  gas  properties  can  be.  studied 
in  the  known  state  of  the  ho*  gas  in  the  shock  tube,  and  consequently 
applied  to  the  flight  situation 

Based  on  this  philosophy,  the  staff  of  the  AVCO  Research 
Laboratory  have  investigated  many  of  the  critical  questions  of  hyper¬ 
sonic  flight  and,  through  exploitation  of  the  straight  shock  tube, 
have  obtained  engineering  answers  which  are  being  used  for  design 
and  development  of  hypersonic  vehicles  ’  The  research  we  have 
conducted  has  been  aimed  at  answering  the  most  critical  questions, 
i.e.  ;  (11  Can  we  learn  enough  about  the  chemical  kinetics  of  the  high 
temperature  air  and  component  gases  to  define  the  thermodynamic 
state  in  both  the  shock  tube  and  flight;  12)  Is  the  aerodynamic  heat 
transfer,  as  determined  from  a  laminar  boundary  layer  theory  in¬ 
cluding  diffusion  of  atoms  a  good  approximation  to  the  high  temper¬ 
ature  situation  or  are  there  unanticipated  phenomena  present;  (3) 

What  IS  the  radiative  emissivity  of  high  temperature  air;  and  (4)  Do 
the  other  high  temperature  gas  properties,  such  as  the  electrical 
conductivity,  complicate  hypersonic  flight? 

Tlie  results  of  these  investigations  and  the  techniques  and 
instrumentation  employed  will  now  be  summarized.  Tltese  descrip¬ 
tions  demonstrate  the.  versatility  of  the  shock  tube  as  a  facility  for 
conducting  hypersonic  research  better  than  any  broad  review  of 
shock  tube  technology  and  speculation  about  its  possible  applications. 


GENERAL  DESCRTFTION’ 


The  schematic  diagram  sho\^’n  in  Fig.  gives  the  notatiori 

and  nomenclature*  of  this  paper  The  shock  tubes  used  in  these  ex¬ 
periments  Lave  been  described  in  »he  literature 

In  general,  they  consist  of  high  pressure  dri'  crs  depending  on  com¬ 
bustion  of  oxygen,  hydrogen  and  h'lium  to  create  the  strong  shock 
waves'*  Systematic  experiments  ha'^e  shown  that  ospecially  for  large 
shock  tubes  the  mixing  of  the  combustible  gases  is  extremely  impor¬ 
tant  and  difficult  If  th»  gases  arc-  properly  mixed  then  controlled, 
non -detonating,  complete  burning  has  been  found  possible  for  mixtures 
of  oxygen  and  hvulrogen  (in  stoichiometric  ra*iol  with  a  diluent  helium 
forming  between  70  and  80%  of  the  total  mass  of  the  mixture  The 
strongGSt  shock  waves  which  can  be  created  by  combustion  of  oxygen 
and  hydrogen  arc  produced  by  mixtures  oopS.aining  (30%  helium,  vvhich 
bum  quite  violently  with  unpredictable  detonations.  Not  much  shock 
strength  is  lost  by  diluting  the  mixture  to  70%  helium 

The  performance  of  constan'  volume  combustion  driv<;rs  is 
quite  predictable  Th-^  speed  of  sound  iu  'b‘=>  hot  dti'  ev  gas  has  been 

measured  from  the  lime  interval  between  the  passage  of  the  expansion 
wave  at  two  pressttvo  gagi?s  m  the,  dri'er  TV.e  rneasuremen’  is  in 
good  agreement  wifh  sound  speeds  calc  ulated  consichring  dissociation 
of  water  vapor  in  the  combustion  products 

Th*  resultant  shock  stiength^,  at:ainable  are  shown  m  Fig  (2) 
The  experimentally  achieved  velocmes  agree  reasonably  with  the 
calculated  predictions  .IIow<ve/  this  is  not  a  jiar*  icul  arly  impor¬ 
tant  feature  because  in  ceneral  the  shock  Mach  numbc?r  is  not  a 
critical  variable  as  it  is  usually  not  required  precisely  Repro¬ 
ducibility  of  a  half  a  Mach  number  is  aieragr  tor  experiments  where 
no  special  emphasis  has  been  placed  on  controlling  this  scatter 

The  speed  of  the  shock  wave  j«  morjhtored  as  i?  travels  down 
the  tube  Tina  mcasuromont  has  been  accomplished  in  a  number  af 
wavs  Shock  waves  are  easily  detected  by  pressure,  heat  transfer, 
gas  conductivity,  denr.ity  and  a  number  of  optical  measuretnents  In 


*  Of  the  many'  different  methods  of  creatine  strong  shock  wave*, 
such  as  healed  light  gases,  double  diaphragms,  def onations,  solid 
explosives  and  electrical  discharge,  unlv  controllable  combustion  ap¬ 
pears  to  be  at  present  at  th"  state  of  development  where  it  can  be  used 
reliably . 
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Fig.  1  —  (a)  Schematic  diagram  o£  shock  tube. 

(b)  X-T  diagram  showing  the  progress  of  the  shock  wave 
and  the  expansion  wave  following  the  diaphragm  burst. 
The  gases  which  were  originally  separated  by  the  dia¬ 
phragm  are  separated  by  the  contact  discontinuity 

(interface). 

(c)  Pressure  distribution  at  several  typical  times  in  the 
shock  tube.  (Time  t^  and  from  (b)). 


Fig.  2  —  Relation  between  shock  Mach  number  and  disphragm  pressure 
ratio  for  three  typical  driver  gases,  (a)  combustion  of  O2, 

and  N2,  (b)  cold  hydrogen,  and  (c)  cold  helium  (Equilibrium 

air  is  the  driven  gas).  Experimentally  measured  values  are 
shown. 
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air,  in  the  r3.age  o£  shock  wave  speetf#  B  to  times  the  speed  of 
sound  at  ro6m  tempeyature,  the  ioiiizaiion- sensitive  picLu.p  which 
measures  the  change  of  resistance  of  an  air  gap,  is  simple  and  ex¬ 
tremely  useful.  Almost  any  degree  of  accuracy  can  be  achieved  by 
properly  choosing  the  time  reeolution  and  pickup  spacing.  A  typical 
oscillograph  is  shown  in  Fig.  (3)  The  sweep  of  a  norrnal  oscilloscope 
has  been  extended  as  a  raster  so  that  very  long  times,  resolvable  to 
very  small  intervals,  are  displayed  Besides  ionization* sensitive 
pickups,  both  heat  transfer  gages  and  phototubes  have  been  used  to 
measure  the  ahock  velocity  in  the  experiments  described 

In  order  to  perform  experiments  in  air  free  of  impuritiea,  it 
has  been  found  advisable  to  flov»  gas  through  the,  low  pressure  shock 
tube  section,  instead  of  using  a  static  test  gas  This  procedure  is 
especially  important  in  spectroscopic  work  It  has  been  found  that  a 
large  amount  of  extraneous  material  can  diffuse  out  of  the  shock  tube 
walls,  "o'‘-rings,,  etc.,  if  a  static,  low  pressure  gas  is  used  A  typi¬ 
cal  flow  system  consists  of  a  throttling  valve,  and  a  check  valve  (which 
is  closed  by  the  pressure  behind  the  shock  wave)  near  the  diaphragm 
in  the  gas  inlet  line,  and  another  check  valve  cift  of  the  station  where 
measurements  are  made  in  a  line  from  a  vacuum  pump  The  flow 
volocity  of  the  test  gas  is  kept  insignificantly  small  compared  to  the 
shock  wave  speed. 

The  properties  behind  a  normal  shock  can  be  calcniated  from 
the  conservation  laws.  i.  e.  ,  the  Rankine  -  Hugoniot  equations  For 
high  temperature  applications  thermodynamic  equilibrium  properties 
of  the  gas,  derived  from  statistical  mechanics,  must  be  substituted 
for  the  Qcmatton  of  s'ate.  For  air.  Hilsenratb  and  Beckett' 

Gilmore^^®'  and  Logan  and  Tre.anor^^^^  have  calculated  the  proper¬ 
ties  of  air  in  thermodynamic  equilibrium.  These  data  have  been  used 
by  Feldman^'^J  to  sob'e  the  normal  shock  and  other  typica.l  shock  tube 
and  hvpersoijic  flight  si'uatioas. 
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Fig.  3  Typical  oscillograph  from  the  shock  velocity  timing  circuit. 

The  oscilloscope  sweep  starts  at  the  top  and  makes  a  complete 
cycle,  left  to  right  and  back,  every  100  microseconds.  The 
bright  spots  are  produced  by  the  output  from  a  crystal  oscil¬ 
lator,  spaced  exactly  10  microseconds  apart.  The  time  of 
arrival  of  the  shock  wave  at  each  successive  ionization-sensi¬ 
tive  pickup  is  denoted  by  the  sudden  break  in  the  sweep  (change 
in  gap  resistance  caused  by  the  rapid  thermal  ionization  of  the 
gas).  Sample  oscillograph  has  signals  from  12  pickups,  spaced 
approximately  a  foot  apart. 
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THERMODYNAMtC  EQUJLLlBJglUM 


Almost  all  shock  tube  meas'Urftfiaeut#  in  high  temperature 
gases  by  nature  of  their  short  steady  State  depend  to  some  degree 
on  the  assumption  that  thermochemical  equilibrium  exists  in  all  of 
the  hot  gas  region.  Much  research  has  been  pointed  at  ascertaining 
this  fact.  This  work  has  mostly  been  directed  toward  measuring 
vibrational  relaxation  time  or  toward  studying  the  dissociation  or  re¬ 
combination  process.  At  the  lowest  temperatures  in  the  region  cov¬ 
ered  by  these  investigations.  Ms  =  b.  0,  the  vibrational  relaxation 
times  in  air  are  already  very  fast.  These  rates  have  been  measured 
by  Blackman^^^)  and  others,  (24)  in  oxygen  and  nitrogen,  and  are 
shown  in  Fig.  (4)  applied  to  air  under  the  assumption  that  the  rates 
applicable  are  those  at  the  total  pressure  of  the  gas,  not  at  the  par¬ 
tial  pressure  of  each  component  It  can  be  seen  that  using  the 
rates  for  the  component  gases,  the  vibrational  relaxation  process  is 
fast  enough  to  leave  the  test  gas  in  equilibrium.  The  data  available 
for  the  relaxation  time  of  either  the  dissociation  or  recombination 
process  ape  not  as  -well  defined  at  the  present  time. 

At  the  AVCO  Research  Laboratory  two  different  types  of 
measurements  have  been  made  to  establish  the  rates  involved  in  the 
dissociation  and  recombination  process.  First,  the  Mach  angles  of 
very  weak  waves  created  by  very  small  surface  irregularities  have 
been  measured  by  Feldman'^'  on  schlieren  photographs.  The  Mach 
number  of  the  flow  is  a  strong  function  of  the  gas  state.  This  meas¬ 
urement  has  been  employed  in  a  number  of  ways.  In  the  original 
experiments  the  Mach  number  in  the  flow  behind  the  moving  normal 
shock  was  measured  by  inserting  a  flat  plate  with  a  sharp  leading 
edge  into  the  flow.  The  experimental  setup  and  representative 
schlieren  photographs  are  shown  in  Fig.  (5).  Results  from  this  ex¬ 
periment  are  plotted  in  Fig.  (6)  .  It  can  be  seen  that  the  Mach 
number  predicted  from  the  equilibrium  assumption  is  actually  meas - 
ured  within  the  scatter  of  the  data.  The  results  from  this  experi¬ 
ment  are  conclusive  but  not  very  quantitative  for  obtaining  rates. 

I 

In  an  attempt  to  measure  the  rate  constant  for  recombination 
of  dissociated  air  precisely,  a  shock  tube  experiment  was  conceived 
to  give  a  very  rapid  expansion  of  a  dissociated  gas.  A  two-dimen-  ( 

stonal  Prandtl- Meyer  expansion,  after  a  compression  on  the  leading 
edge  of  a  wedge,  was  used.  in  order  to  gain  the  maximum  change 
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Fig.  4  —Relaxation  times  behind  a  moving  normal  shock  wave  in  a 
shock  tube.  The  reiaxation  times  (for  vibrational  relaxation 
of  O2  and  from  Ref.  (2b)  and  a  dissociative  relaxation 
baaed  on  a  recombination  rate  constant,  ky,  of  10^^  cm^  sec”^ 
mole"^  and  the  equilibrium  properties  of  air  )  are  given 
in  terms  of  the  length  of  shock  tube  required  for  a  particle 
to  be  immersed  in  the  hot  flow,  region  2  in  Fig.  (1),  for  a 
loiigth  oi  time  equal  to  the  eelaaeticia  ttnris. 


10  ^  a  lb  «e 

SHOCK  MACH  NUMBER- M, 


MODEL  SCRIBED  WITH  HEAVY  LINES  IDENTICAL  TO  EACH  OTHER 


P,  =  I  CM  OF  mercury;  T,  =300°  K;  Mj  =13.85 


I 

I 
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(b) 


rig.  5  —  (a)  Experimental  arrangement  of  model  in  flchlieren  test 
section  i-ti  ehock  tube  diemical  kinetic  ei!;peTijments. 

(b)  ScLlieren  photograph  of  flow  over  modal  showing  Mach 
lines  and  direction  of  line#  predicted  from  equilibrium 
as  sumption. 
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FLOW  MACH  ANGLE  (DEGREES) 
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Fig.  6  •—  Estimates  of  flow  Mach,  number  downstream  of  moving 

shock.  (Relaxation  phenomena  are  includad  when  computing 
sound  velocities).  Gompaxison  of  experiments  with  frozen 
flow  (  /  s  7/5),  equilibrium  flow,  and  the  Mach  number 
which  the  flow  would  have  if  vibration  were  in  equilibrium 
but  no  dissociation  had  occurred. 
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of  Mach  number  between  the  extreme  chemical  processes,  i. , 
equilibrium  and  "frozen"^  flow,  the  ivedge  angle  was  adjusted  to  give 
very  nearly  sonic  flow  on  itp  forward  edge.  Various  turning  angles 
of  the  expansion  were  used.  The  general  concept  of  the  experiment 
is  that  if  recombination  of  atome  takes  a  long  enough  time,  then 
small  disturbances  aft  of  the  cornar  will  create  Mach  waves  which 
index  the  local  state  of  the  gas.  Unfortunatelv,  all  measurements  of 
the  Mach  angle  even  from  the  first  disturbance  behind  the  corner, 
have  agreed  with  values  which  are  predicted  from  the  equilibrium 
assumption,  as  shown  in  Fig.  (7l.  As  a  result,  a  rough  estimate  of 
a  lower  limit  for  the  recombination  rate  constant  is  the  only  informa¬ 
tion  that  has  been  obtained  from  this  experiment  to  date. 

As  a  by-product  of  this  experiment  it  was  observed  that  the 
oblique  shocks  on  the  forward  portion  of  the  wedge  assumed  an  angle 
predictable  from  equilibrium  theory.  This  condition  is  very  nearly 
the  same  as  at  the  stagnation  point  of  a  blunt  body.  Th'^  comparison 
of  the  observed  and  calculated  oblique  shock  angles  are  shown  on 
Table  (1*.  Details  of  this  work  have  been  reported  in  Ref,  (9). 

The  second  method  used  by  Camm  and  Keck^^^'  to  observe  the 
rate  processes  i«  a  system  which  monitors  the  radiative  relaxation 
behind  a  moving  normal  shock.  The  radiation  behind  a  moving  normal 

shock  IS  observed  through  a  vertical  slit  collimated  'o  produce  time 
resolution  better  tlian  0„  5  /^sec.  ,  when  the  signal  is  observed  by  a 
photomultiplier.  The  radiation  has  also  been  dispersed  by  a  mono.- 
ebroruator  and  viewed  over  several  limited  wave  length  bands. 

Results  have  been  obtained  in  oxygen  and  nitrogen.  The  re¬ 
combination  rate  constant  kr,  at  or  near  equilibrium,  obtained  from 
these  experiments  appears  to  be  in  general  agreement  with  the  re¬ 
sults  of  Glick  and  Wurster(24l^  ^  typical  oscillograph  showing  the 
radiation  relaxation  is  presented  in  Fig,  (8).  The  rate  constant,  kr, 
(in  cni^  sec.”^  moles"^'  estimated  to  be  of  the  order  of  10^^  for 
the  component  gases  and  essentially  independent  of  temperature  by 
several  investigators,  appears  to  be  approximately  as  estimated. 


*  The  "frozen**  and  "equilibrium"  situation  can  best  be  distinguished 
by  considering  the  relaxation  time  for  recombination  of  atoms  Z~m  as 
compared  to  a  characteristic  flow  time,  TV  •  The  frozen  flow  is  de¬ 
fined  as  ‘TV  and  equilibrium  as  »  Tm. 


Fig.  7  —  Comparison  oi  experimental  Mach  lines  with  equilibrium 
eind  frozen  Mach  angles,  (l^equil  =  equil  = 
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-H  H-|;iSEC 


BLUE 


Pbotomultlpliar  tr*c*a  from  xadiativa  relaxation  «xperimsat. 

Tha  rfuUatloa  iA  the  red  ead  blue  bande  Indicatee  the  eeine 
xelexetion  time  even  though  one  contaius  CN  radiation  end 
the  other  not.  Luxninoue  front*  have  previouely  been  attributed 
\o  CN  redixtioa  la  other  gaeoe. 


The  relaxation  rates  iti  ai*  are  likely  to  be  considerably  faster 
than  they  are  for  the  component  gases,  due  to  the  possibility  that 
recombination  is  not  a  simple  three-body  process,  but  proceeds  by 
intermediate  steps.  As  a  result  definitive  results  will  not  be  avail¬ 
able  until  direct  measuremente  have  been  made  in  air.  Indications* 
are  that  the  rates  are  at  least  a  factor  of  ten  better  than  the  10^^ 
measured  for  the  component  gases.  A»  a  reeult  a  recombination  rate 
constant  of  10 has  been  used  in  a  calculation  to  determine  the  dis¬ 
sociation  relaxation  time  in  the  shock  tube.  The  results  are  shown 
in  Fig.  (4)  and  imply  that  the  gas  in  the  shock  tube  is  in  chemical 
equilibrium. 


*  Some  of  these  indications  are,  besides  the  direct  experiments 
described  above,  heat  transfer,  electrical  conductivity  and 
radiati\c  cmissiv-ity  measurements,  all  independent  asid  abso¬ 
lute,  which  each  agree  only  with  theoretical  predictions  based 
on  the  fact  that  the  air  behind  the  moving  nornnal  shock  is  in 
equilibrium 
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f HYPERSONIC  AERODYNAMICS 


The  inviscid.  flo^iv  At  hypersonic  7>iacli  nuinbers  is  one  problem 
for  ^bich  the  shock  tube  is  severely  limited.  Still  there  are  advantages 
because  in  the  flow  about  blunt  bodies,  where  eesentially  all  the  air 
which  enters  the  boundary  layer  went  through,  a  strong  normal  shock, 
the  inviscid  flow  Mach  number  does  not  exceed  valuee  around  three 
even  at  Mach  20.  As  a  result  the  shock  tube  boundary  layer  ie  never 
radically  different  from  flight.  However,  the  straight  shock  tube 
cannot  give  exactly  correct  answers  for  pressure  distribution  at  the 
proper  flow  Mach  number  or  density  ratio  across  the  bow  shock. 

Out  effort  in  this  field  has  been  limited  because  other  types  of  facil¬ 
ities,  such  as  ballistic  ranges,  by  the  use  of  optical  techniques  in¬ 
volving  measurement  of  Mach  number  distribution  from  Mach  lines, 
appear  to  be  suited  for  producing  significant  results.  In  the  shock 
tube  it  will  suffice  to  measure  the  pressure  distribution  to  define  the 
inviscid  flow  parameters  which  exist  during  experiments. 

Boundary  layer  problems  which  can  be  separated  from  the 
inviscid  flow  are  particularly  suitable  for  investigation  in  shock  tubes. 
The  first  and  most  serious  problem  studied  was  the  heat  transfer 
through  a  laminar  boundary  layer  at  the  stagnation  point  of  blunt 
bodies^®)'  (10),  The  stagnation  point  of  a  blunt  body  has  several 
advantages  which  make  it  particularly  attractive  for  heat  transfer 
experiments.  First,  the  most  severe  heating  occurs  at  this  point; 
secondly,  direct  simulation  of  the  flight  conditions  is  possible;  and, 
thirdly,  the  problem  of  conductive,  convective,  and  diffusive  (diffu¬ 
sion  of  atoms  as  differentiated  from  molecules)  heat  transfer  has 
been  treated  theoretically  at  the  laboratory  by  Fay  and  Riddell 
probably  the  only  example  of  a  boundary  layer  solution  which  is 
capable  of  including  non-equilibrium  gas  effects. 

Considering  this  problem,  it  can  be  shown  that  the  heat 
transfer  rate  through  a  laminar  boundary  layer  at  the  stagnation 
point  at  low  speeds,  where  the  product  of  the  density  and  viscosity, 

(  ),  in  the  boundary  layer  does  not  vary  from  the  free  stream 

to  the  wall  and  density  varies  inversely  with  temperature (^  ^),  can 
be  expressed  in  a  non-dimensional  form,  by 

*0763  0^^"*  (1) 
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wile  re  Nu  ,  Se  .  O'  are  the  Nusselt,  Reynolds, 

and  Prandtl  numbers,  respectively. 

The  numerical  aolution  of  the  equilibrium  boundary  layer 


problem  treated  by  Pay  and  Riddell 


(10> 


can  be  correlated  by  con* 


side  ring  the  ratio  of  (  )  at  the  edge  of  the  boundary  layer  to 

( J  wall,  follows; 

^04. 


(2) 


where  is  the  Lewis  number,  and  is  the  fraction  of  the  total 


enthalpy  which  is  in  dissociation. 


From  this  expression  the  heat  transfer  rate, 
stagnation  point  can  be  expressed  by  ® 


at  the 


r 

.  1/ 

i,.  0.763 

^)s  L* 

cr'^'VidxL] 

(3) 


where  is  the  enthalpy  of  the  gas  at  the  wall  tempera- 


ature  and 

point. 


'^/i% 


x*o 


is  the  velocity  gradient  at  the  stagnation 


It  can  be  seen  from  Eq.  (3)  that  simulation  of  heat  trans¬ 
fer  requires  that  the  conditions  at  the  edge  of  the  boundary  layer 
are  properly  reproduced,  that  the  wall  conditions  be  similar,  and 
that  the  velocity  gradient  at  the  stagnation  point  have  the  same 
value  in  a  model  test  as  in  the  flight  situation  simulated.  In  addi¬ 
tion,  the  chemical  state  of  the  air  in  the  boundary  layer  must  be 
considered. 

The  solutions  of  the  Rankine-Hugoniot  equations.  Ref.  (7) 
and  Fig.  (3),  make  it  evident  that  the  stagnation  density  and  en¬ 
thalpy  encountered  in  hypersonic  flight  are  easily  accessible  in 
ordinary  shock  tubes.  In  addition,  material  limits  demand  that  the 
wall  temperatures  be  similar  in  the  two  situations. 

The  velocity  gradient  at  the  stagnation  point  is  derived  from 
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the  pressure  distribution  In  hypersonic  aerodynamic#  the  !^ewtoni4n 
estimate  has  be&n  found  to  predict  the  correct  pressure  distribution 
very  closely  even  at  Surpirisingly  low  Mach  nutnbera  for  blunt  bodies. 
This  estimate  is  independent  of  the  Mach  number  of  the  flow.  Opti¬ 
cal  measurement  of  Mach  angles  in  the  flosy  from  fchlieren  photo¬ 
graphs  have  been  used  mostly’^'  From  the  Mach  number  dist'ibution 
the  pressures  can  be  calculated  by  following  an  isentropic  expansion 
on  a  graphical  representation  of  the  equilibrium  properties  of  air. 
Such  as  the  Mollier  diagram  of  Ref.  (7)  A  typical  Schliercn  photo¬ 
graph  and  the  pressure  distribution  calculated  therefrom  are  shown 
in  Fig  no)  The  pressure.s  are  compared  to  the  Newtonian  estimate 
and  show  reasonable  agreement  Tn  particular  the  pressure  near  the 
sonic  point  is  in  good  agreement  As  a  result  it  is  concluded  that 
the  velocity  gradient  at  the  stagnation  point  can  be  talten  as  the 
N  ewto.iian  value  in  the  shock  tube  as  well  as  in  hypersonic  flight 
and  therefore  is  independent  of  the  Mach  number  in  front  of  the  bow 
shock 


As  a  resull  of  these  arguments,  it  is  seen  that  the  shock 
tube  car  produce  the  environment  encountered  at  the  stagnation  point 
in  hypersonic  flight  The  initial  shock  tube  conditions,  and  the 
flight  velocity  and  altitude  at  which  the  same  stagnation  enthalpy 
and  density  exists,  can  be  specified  h'ig  (11)  is  a  plot  of  such 
Simulation.  The  region  covered  by  the  stagnation  point  measurements 
described  is  outlined  Trajectories  for  a  number  of  contemplated 
bypersonice  vehicles  have  been  added  to  this  chart  from  Ref.  (li). 
The  straight  shock  tube  is  seen  to  cover  much  of  the  critical  heat 
transfer  region 

The  considerations  outlined  so  far  are  rompleto  if  tho  flow 
chemistry  is  similar,  i.  e  ,  if  the  boundary  layer  in  both  the  shock 
tube  and  flight  is  either  frozen  or  in  equilibrium  In  small  scale 
tests,  the  boundary  layer  is  much  thinner  than  it  is  on  the  larger 
body  and  an  atom  therefore  spends  less  time  diffusing  from  the 
boundary  layer  edge  to  the  wall  Aa  a  result  there  are  situations 
in  which  the  chemistry  i«  not  exactly  simulated 

The  case  where  the  boundary  layer  is  neither  in  equilibrium 
nOr  completely  frozen,  but  has  a  finite  recombination  rate  for  atoms, 
has  been  analyzed  by  Fay  and  Riddell  In  this  reference  it  is 

sho  wn  that  the  laminar  boundary  layer  problem  can  be  solved  by 
considering  a  recornbination  rate  parameter  A  constant  value  of 


*  In  our  high  Mach-  number  shock  tubes,  the  testing  time  is  so  small, 
Fig.  (9),  that  the  requirements  on  a  pressure  transducer  are  very 
severe. 
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MICRO- SECOMDS  PER  FOOT 


Fig.  10  —  (a)  Schlieren  photograph  of  the  flov/  over  the  nose  of  a 

hemisphere-cylinder  model  in  a  1-1/2  inch  I.  D.  shock 
tube.  Note  fine  Mach  lines  whose  inclination  to  the 
body  is  accurately  measurable. 


(b)  Pressure  distribution  calculated  from  Mach  line  measure¬ 
ments  compared  to  the  Newtonian  pressure  distribution. 
The  agreement  is  reasonable,  especially  on  the  aft 
portion  of  the  body  and  near  the  sonic  point. 

Stagnation  conditions  are  equivalent  to  a  flight  velocity 
of  approximately  13,000  ft/sec  a.nd  an  altitude  of  25,  000  ft. 
’Tie  flow  Mach  number  is  approximately  2.  3. 
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ALTITUDE  -  FTx  10* 


m 


24C 


-15,000  FRS. 


Fig.  IX  —  Shock  tuhe  conditions,  Ms  «nd  Pi,  compared  to  Oj>e  flight 
velocity  and  altitude  at  which  the  same  atagnation  enthalpy 
and  preasuze  are  encountered.  The  outlined  region  ie 
covered  by  the  preaent  experimenta. 

Typ  leal  trajcctorieo  of  several  types  of  hypersonic  vehicles 
are  ahoiim  from  Ref.  il2). 


-26- 


this  parameter  implies  that  in  addition  to  the  stagnation  enthalpy,  the 
flow  chemistry  in  the  boundary  layer  is  sinnulated.  At  «*.  given  flight 
velocity,  and.  therefore  stagnation  temperature,  the  recombination 
rate  parameter  is  proportional  to  the  square  of  the  density  and  a 
model  body  radius.  Experiments  performed  on  a  scale  model  at  a 
given  density  simulation,  actually  produce  the  density  which  occurs 
in  flight  at  much  hither  .altitudes,  depending  On  the  scaling,  even  if 
the  boundary  layer  not  completely  in  equilibrium.  The  altitudes 
at  which  the  chemistry  is  simulated  for  a  1/1000  scale  model  are 
shown  in  Fig  121 

In  order  to  make  significant  heat  transfer  measurements ,  an 
instrument  capable  of  microsecond  response  was  required  The  most 
promising  schemes  considered,  both  variations  of  the  resistance 
thermometer  principle,  were  HJ  surface  temperature  measurements 

as  used  by  the  Princeton  group'  '  and  others'  ,  and  (Z)  calorim¬ 
etric  measurements  devised  and  introduced  by  this  laboratory. 

From  considerations  of  the  one -dimensional  heat  conduction 
problem  for  a  finite  metal  thickness  backed  by  a  semi-infinito  elec¬ 
trical  insulator,  die  heat  flux  into  the  gage  for  both  these  schemes 

can  be  found 

The  two  extromes  for  the  ratio  of  the  gage  thickness.,  Jt  ,  to 

the  characteristic  diffusion  depth,  >  yield  the  results  for 

the  two  types  of  gages  If  the  gage  thickness,  £  .  is  very  small 

compared  to  Si  in  the  metal  then  the  gage  has  essenti  ally  no  heat 
capacity  and  assumes  the  surface  temperature  of  the  insulator^i^^^^i^^^^^)- 
A  constant  heat  flux,  CJ  ,  gives  a  voltage,  E{tl,  which  is  proportional 
to  the  instantaneous  surface  temperature.  The  shape  of  E(t)  is 
simply  a  pai-abola  since  the  surface  temperature  rises  with  the  square 
root  of  time. 

E±l 

/t  (4) 

In  the  other  extreme,  where  the  thickness,  X  *  i*  much 
larger  than  the  characteristic  diffusion  depth,  the  temperature  at 
/if  :  ^  19  not  affected  during  a  test  so  that  no  heat  is  transferred 
from  -  1,-^%  to  In  this  case  the  total  beat  input  must 

be  reflected  m  the  average  temperature  rise  and  the  heat  capacity  of 
the  gage.  The  gage  integrates  the  total  heat  input,  therefore,  the 


^  const. 
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MODEL  SCALE  l/tOOO 
FULL  SIZE 


P,  =  0  01  C'T)  Hg 


■Hi 


FLIGHT  VELOCITY -FT/SEC  X  JO 

Fig.  12  •—  Shock  tvihe  Mg  and  compared  to  the  Jiypereonic  flight 
cocditioa  at  vhich  atagsiation  point  enthalpy  and  flow 
ehemlatxy  are  identical  for  a  model  scale  of  1/1000. 


voltage  output  varies  directly  proportional  to  the  total  heat  input.  The 
instantaneous  heat  flux  is  determined  from; 


-i  =  const.  Y—  (5) 

dt 

Both  the  above  analyses  can  be  extended  to  the  case  where  the 
heat  flux  varies  with  time,  i.  e.  ,  q  =  q{t). 

Each  type  of  gage  has  limitations  which  must  be  considered.  The 
thin  gage  cannot  be  so  thick  that  it  has  significant  heat  capacity  and 
causes  surface  temperatures  noticeably  different  from  the  one-dimen¬ 
sional  semi-infinite  case.  The  thick  gage  in  turn  must  be  thick  enough 
so  that  the  heat  loss  from  the  rear  of  the  gage  is  small.  The  detailed 
analysis  and  investigations  made  to  evaluate  the  two  heat  transfer  gage 
principles  are  given  in  Ref  (13) 

The  advantages  of  the  calorimeter  gage  are  that  much  higher 
total  heat  inputs  can  be  handled  without  creating  excessive  surface 
temperatures  and  that  the  data  reduction,  especially  for  time  depend¬ 
ent  heat  transfer  rates,  is  much  simpler.  In  addition,  the  calorime¬ 
ter  gage  requires  no  calibration  because  the  bulk  properties  of  chemic¬ 
ally  pure  metals,  such  as  platinum,  are  well-known.  For  the  thin 
gage  a  calibration  technique  in  which  a  known  electrical  energy  pulse 
of  the  same  time  duration  and  energy  content  is  used  to  determine  the 
constant  in  Eq.  (4) 

Both  types  of  gages  have  been  made  as  simple  as  possible 
because  at  the  stagnation  point  of  blunt  bodies  very  strong  shock  waves 
produce  severe  erosion,  and  destroy  the  surface  and  gage  element. 
However,  the  damage  occurs  only  after  the  contact  surface  has  passed, 
as  seen  from  the  smooth  variation  of  resistance  with  time,  for  both 
the  thin  and  the  calorimeter  gages.  Typical  oscillographs  for  both 
types  of  heat  transfer  gages  are  shown  in  Fig.  (13). 

The  results  of  the  stagnation  point  heat  transfer  experiment 
are  shown  in  Fig.  (14).  The  experiment  was  conducted  at  three  con¬ 
venient  initial  pressures  in  the  low  pressure  section,  pi  m  10,  1.  O, 
and  0.  1  cm  of  Hg.  ,  corresponding  to  flight  altitudes  of  approximately 
20,000,  75,000,  and  120,000  feet.  The  velocity  range  covered  by 
these  experiments  extends  from  velocities  at  which  essentially  no 
dissociation  occurs  to  velocities  beyond  the  satellite  speed,  about 
26,000  ft/ sec.  ,  where  more  than  60%  of  the  energy  is  inv’ested  in 
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Fig.  13  —  Typical  oscillographs  and  data. 

On  all  three  oscillographs  time  increases  from  left  to 
right.  The  sweep  speed  is  10  /x  sec/interv'al.  The  oscillo^ 
scope  is  triggered  before  the  moving  shock  arrives  at  the 
stagnation  point.  The  break  in  the  horizontal  trace  denotes 
the  arrival  of  the  shock  at  the  gage.  The  break  of  the 
curve  accompanied  by  the  end  of  the  constant  heat  transfer 
denotes  the  arrival  of  the  interface. 


(a)  Record  from  a  thin  heat  transfer  gage  at  a  relatively 
low  shock  Mach  number  (Note  parabolic  signal  rise) 
fb)  Record  from  a  thin  heat  transfer  gage  at  a  higher  shock 
Mach  number.  Note  that  the  extreme  conditions  at  the 
high  stagnation  enthalpy  degrade  the  quality  of  the  signal, 
(c)  Calorimeter  gage  record.  (Note  linear  signal  rise). 
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diB.«ocl«tior)  energy  eit  the  et«ign»tion  point.  The  theofetical  curves 
are  plotted  izom  the  equation  which  correlates  the  equilibrium  boundary 
layer  solutions  of  Fay  and  Riddell  Sq.  (3).  The  L.ewis  number 

has  Wen  estimatod  and  taken  to  be  1.4  and  the  Prandtl  number  0.  71. 
The  theory  and  data  can  be  seen  to  agree  within  the  scatter  of  the 
data.  An  uncertainty  in  the  absolute  value  of  the  theoretical  predic¬ 
tion,  due  to  the  uncertainties  in  predicting  the  viscosity  and  Lewis 
number  of  the  hot  gas.  is  of  the  same  order  as  the  scatter  of  the  data. 

As  discussed  previously  the  significcutce  of  this  data  can  be 
extended  to  higher  iiltitudes  if  simulation  of  the  boundary  layer  chem¬ 
istry,  based  on  the  recombination  parameter  is  considered.  The  best 
available  estimates  of  the  recombination  rate,  such  as  the  data  of 
Ref.  (9)  and  (23),  introduced  into  the  non-equilibrium  boundary  layer 
solution  of  Fay  and  Riddell  indicate  that  the  boundary  layer  is 

neither  in  equilibrium  nor  completely  frozen  for  the  experiments  per¬ 
formed  at  initial  pressures  of  0.1  cm  of  Hg.  However,  the  total  heat 
transfer  by  conduction  and  convection  is  very  nearly  the  same  for  the 
frozen  and  the  equilibrium  boundary  (at  Lewis  number  of  1.  4)  except 
in  the  extreme  case  where  the  boundary  is  completely  frozen  and  the 
wall  is  non-catalytic  to  the  atom  recombination.  Because  for  present 
recombination  rate  estimates,  the  boundary  layer  is  still  only  partially 
frozen  at  *  0.1  cm.  ,  the  present  experiments  are  not  capable  of 
defining  the  flow  chemistry  from  heat  transfer  measurements.  These 
effects  are  discussed  in  detail  in  Ref  flO). 

If  the  air  were  not  to  equilibrate  across  the  bow  shock  at  the 
Stagnation  point,  the  stagnation  pressure  would  be  reduced  and  the 
heat  transfer  lowered  correspondingly.  At  the  extreme  of  the  condi¬ 
tions  covered  by  these  experiments,  this  would  amount  to  a  reduction 
of  20%  in  the  heat  transfer  rate.  However,  in  view  of  the  evidence 
presented  earlier  indicating  that  equilibrium  is  attained  on  the  leading 
edge  of  a  two-dimensional  wedge,  the  gas  state  is  not  eericuely  in 
doubt  Further,  it  was  shown  that  thermodynamic  equilibrium  is 
reached  across  a  moving  shock  and  therefore  equilibration  across  the 
bow  shock,  where  e\'en  higher  temperatures  and  preBsurea  exist, 
will  be  much  more  rapid 
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EMISSIVITT  of  HIGK  VESIPSRATURE  AIR 


A  second  major  investigation  (Undertaken  by  Keck,  KiveS  and 
Wentink^^'^),  was  the  determination  of  the  radiative  emissivity  of  high 
temperature  air.  Prior  to  the  performance  of  these  experirnents, 
rough  estimates  were  made  by  H.  Bethe  and  H.  Meyer  at  AVCO  and 
others  elsewhere  which  indicated  that  the  NO  ^  and  (f  bands  might 
have  a  transition  probability  (f  number  of  the  order  of  0.1)  which 
would  give  a  considerable  amount  of  radiant  heat  transfer.  In  fact, 
radiative  heat  transfer  was  estimated  to  be  roughly  of  the.  same  order 
as  the  aerodynamic  stagnation  point  heat  transfer  for  a  one  foot  nose 
radius  sphere  flying  at  Mach  20  at  100,  000  ft  altitude. 

Because  radiation,  for  optically  thin  layers  with  small  emis- 
sivities,  scales  in  direct  proportion  with. the  thickness  of  the  radiating 
gas  layer,  the  direct  model  test  technique  is  not  applicable,  mainly, 
because  the  radiative  energy  would  be  extremely  low  and  difficult  to 
detect.  It  was  decided  to  study  the  radiation  from  a  layer  of  gas 
shock-heated  to  the  temperature,  density  and  thermodynamic  state  of 
interest  to  hypersonic  flight  problems  by  reflecting  a  strong  shock 
wave  from  the  closed  end  of  a  shock  tube.  The  extension  to  flight 
configurations  remains  a  geometric  problem.  A  quartz  window  was 
inserted  in  the  end  of  the  lube  for  observation.  The  original  meas* 
urements  were  performed  by  using  the  "thin"  heat  transfer  gage, 
as  used  for  aGrodynamic  heat  transfer  measurements,  behind  a  quartz 
window  These  experiments  established  the  order  of  magnitude  and 
temperature  dependence  of  the  integrated  radiative  energy  flux. 

For  identification  of  the  radiating  species  and  quantitative 
emissivity  measurements,  time  resolved  spectrograms  and  calibra¬ 
ted  monochromator-photocell  measurements  were  made.  The  spec¬ 
troscopic  work  was  devoted  mainly  to  identification  and  measurement 
of  relative  intensities  of  the  radiating  species  Time  resolved  emis¬ 
sion  spectra  were  obtained  in  order  to  separate  the  radiation  which 
occurs  eifter  the  driver  gases  (or  mixed  region)  encounters  the  re¬ 
flected  shock.  Time  resolution  of  about  10  yta  sec.  has  been  used. 

The  expe r iTnenta.1  setup  is  shown  in  Fig.  {15)  and  a  typical 
spectrum  in  the  region  of  2700  to  4000  ^  is  shown  in  Fig.  (1^)- 
The  second  positive  b«nd  system  of  nitrogen  and  CN  are  clearly 
detectable.  Especially  at  wave  length*  smaller  than  3000  the 
NO  band  systems  are  expected  to  have  appreciable  inten  »ity 
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Fig.  i5  Ebq>eriinftntal  arrangement  for  radiative  emieeivity  measure 
mezits,  showing  end  of  ehock  tube,  window,  epectrograph 
end  drum-  camera* 
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Fi«.  16  —  Time  resolved  spectra  in  the  region  of  2700  -  4200  A.  Note 
the  build-up  oi  radiative  intenaity  as  the  gas  layer  behind  the 
reflected  ahock  grows,  and  the  end  of  the  "clean"  radiation  at 
approximately  30  /x  sec  when  the  driver  gas  and  reflected 
ahock  meet. 


Positive  idftfttification  of  NO  bands  at  SOOO^Kis  difficult  because  of 
the  lack  of  characteristic  structure, 

in  the  visible  and  near-infsrarftd  ragion,  the  radiation  shows 
little  etructure.  The  radiation  above  5300  A  can  best  be  attributed 
to  the  First  positive  system  of  nitrogen.  The  CN  violet  and  red 
band  systems  are  also  visible.  The  lack  of  structure  in  the  whole 
region  may  be  due  to  the  reaction  of  O  4-  electron—^  O  h>/,  however 
the  continuous  distribution  makes  identification  difficult. 

Absolute  measurements  of  the  spectrad  emissivity  of  high 
temperature  air  in  the  region  of  Z300  ^  to  9500  have  been  nnade, 
using  a  calibrated  monochromator  and  photocell.  The  measurements 
were  made  at  a  shock  Mach  number  of  13.2  t  0.  5  and  an  initial  pres¬ 
sure  of  1  cm  of  Hg .  ,  resulting  in  an  equilibrium  temperature  of 
8050**K  and  a  density  of  .  85  i  .  04  normal  atmospheric,  sea  level 
density.  The  measured  wave  length  dependence  of  radiation  from  air 
at  a  temperature  of  8050°K  is  plotted  in  Fig.  (17).  The  significant 
contributing  species  are  tabulated  in  Table  II.  where  the  f  numbers 
and  the  emissivity  of  a  one  cm.  thick  layer  are  given.  In  this  table 
positive  spectroscopic  identification  can  be  claimed  for  the  CN  and 
N2  band  systems. 


TABLE  II 


Emissivity  of  optically  thin  air  attributed  to  six  radiating  species. 

The  f  numbers  except  for  CN  are  deduced  from  absolute  emissivities 
at  appropriate  wave  lengths.  They  are  subject  to  the  assumption 

that  the  large  number  of  lines  in  the  molecular  bands  can  be  treated 
as  a  continuum.  They  are  defined  as  the  f  numbers  per  initial  elec¬ 
tronic  State.  (Air  at 8000°K  and  //^  ::  0.85) 

Source 

f 

number 

%  emissivity  cm~^ 

e/L 

NO  ( >5  and  r  )  (.  2  /U  c  A  ) 

.  0025 

1. 1 

N2  (2  nd.  pos.  } 

.  07 

1.  i 

CN  (violet) 

.  1 

.  1 

0~  (free-bound) 

.  4 

N2  (1st.  pos  )(A<  lyW) 

.  02 

.  6 

0  ♦  ft  (free-free)  (  X  <  1 /U  ) 

.  1 

Total  Theoretical  Ernissivity  ( 

2  <A  <  1/; 

) 

3.  4 

Total  Experimental  Emissivity 

(.2  <  A<l/« 

3.5 

SPECTRAL  STERRAOIANCY  W/ 


Fig.  17  —  Monochromator-photocell  measuremont®  of  the  distribution 
and  absolute  intensity  of  the  radiation  from  1  cm  of  air 
at  approximately  SOOO°K  and  .  85  compared  to  the 

estimated  radiation 
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For  the  continuum  radiation,  the  theoretical  conwiderations  oi  Kivel 
coupled  with  the  general  fe^atures  of  the  spectra  niake  these  emis- 
sivities  plausible*.  The  theoretical  analysis  (^7)^  result*  of  which  are 
shown  in  Fig-  (17^  together  with  the  meaaurennent*,  has  as  free  par¬ 
ameters  only  the  f  number  of  the  bands  of  N2  anti  NO.  The*e  have 
been  adjusted  to  fit  the  spectral  distributions  measured  in  the  appro¬ 
priate  region  of  the  Spectrum. 

It  is  important  to  notice  that  the  f  number  for  NO  is  consid¬ 
erably  less  than  the  estimated  0.  i.  In  addition,  the  total  emissivity 
is  approximately  .  050/cm  at  8000®K  and  0.  85.  Measurements 

of  the  temperature  and  density  dependence  of  the  radiation  were  made 
using  the  monochromator  and  photocells  The  resultant  temperature 
and  density  dependence  of  the  absolute  emissivity  per  cm.  of  air  is 
shown  in  Fig.  (18).  The  results  of  Fig.  (18)  can  he  used  to  compare 
the  heating  at  the  stagnation  point  due  to  radiative  heat  transfer  to 
the  laminar  aerodynamic  heat  transfer  rates.  It  can  be  shown  that 
in  flight  at  a  Mach  number'  of  20  at  approximately  120,  000  ft.  ,  the 
radiative  heat  transfer  for  a  one  foot  radius  sphere  will  be  only  about 
10%  of  the  aerodynamic  stagnation  point  heat  transfer 


*  The  CN  and  free-free  and  free-bound  radiation  associated  with 
Oxygen  have  been  estimated  based  on  other  data. 
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DENSITY  RELATIVE  TO  NORMAL 


Fig<  18  — •  Temperature  and  density  dependence  of  tLe  emissivity  of 

&ir  from  the  theoretical  calculatlone  of  Ref.  (27).  These 
eetimatee  have  been  checked  in  a  limited  region  by  the 
QiqpeTimente  deecribed* 
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ELECTRICAL  CONDUCTIVITY  OF  HIGH  TEMPERATURE  AIR 


Another  investigation  carried  on  at  this  laboratory  i» 
work  of  S.  G.  Lin  a  study  of  the  electrical  conductivity  of 

thermally  ioniced  air. 

The  techniquea  employed  in  this  experiment  are  similar  to 
those  reported  in  detail  in  Ref.  H)*  An  axisymmetric  magnetic  field 
ie  produced  in  a  glass  shock  tube  section.  When  the  shock  heated 
gas.  which  is  made  conducting  due  to  thermal  ionization,  arrives  at, 
and  enters  the  stationary  magnetic  field,  the  magnetic  lines  of  force 
are  displaced.  By  measuring  this  displacement  of  the  magnetic  field 
through  the  use  of  a  small  pick-up  coil  nearby,  the  electrical  con¬ 
ductivity  of  the  gas  can  be  deduced  from  the  known  gas  velocity  and 
the  strength,  geometry  and  displacement  of  the  magnetic  field  It 
is  convenient  to  calibrate  the  system  by  shooting  a  slug  of  metal  of 
known  conductivity  through  the  field  at  a  known  speed.  A  schematic 
representation  of  the  experimental  arrangement  is  shown  in  Fig.  (19). 

Typical  oscillogram  traces  from  the  pick-up  coil  signal  of 
this  experiment  performed  at  a  shock  Mach  number  of  12.  0  are  shown 
in  Fig.  (201,  along  with  the  signals  obtained  in  the  calibration  using 
an  aluminum  bar.  From  the  amplitude  of  the  voltage  peak  and  the 
velocity  of  both  calibration  and  experiment,  and  also  the  known  con¬ 
ductivity  of  aluminum,  the  gas  conductivity  can  be  calculated. 

It  can  be  seen  from  the  oscillogram  traces  that  the  positive 
and  negative  pulses  are  the  same  amplitude  in  the  shock  tube,  as 
well  as  for  the  aluminum  bar,  which  implies  that  the  electrical  con¬ 
ductivity  distribution  must  be  quite  uniform  throughout  the  hot  gas 
region.  Furthermore,  since  the  experimental  pulse  shape,  normal¬ 
ized  to  the  same  sweep  speed  as  the  calibration,  is  very  similar  to 
the  aluminum  result,  the  electrical  conductivity  and  hence  the  degree 
of  ionization  must  have  built  up  very  rapidly  behind  the  shock  front. 
Since  the  electrical  conductivity  or  degree  of  ionization  behind  the 
ehock  front  built  up  very  sharply  at  all  shock  Mach  numbers  inves¬ 
tigated  (M*  r  10  to  18),  it  may  be  concluded  that  the  equilibrium  de¬ 
gree  of  ionization  was  reached  aimoet  immediately  and  maintained 
for  all  these  experiments 

The  measured  electrical  conductivity  in  practical  units  ( Mhos  /  cm 
is  plotted  in  Fig.  (21)  compared  to  theoretically  calculated  conductivities. 
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Fig.  19  —  Sch.«maUc  di*grAx»  •hewing  the  principal  elements  of  the 

shock  wave^magnetic  field  interaction  experiment. 
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Fig.  20  — »  Comparison  o£  typical  oscillograms  from  the  shock 
wave  and  from  the  aluminum  bar  calibration  In  the 
magnetic  field  interaction  experiment. 
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SHOCK  MACH  HQ, 


Fig.  21  Comparison,  oi  xnonsured  electrical  conductivity  of 
sh  ock-hcnted  nir  vlth  theory  baeed  on  et^uiiibrium 
degree  of  ionization,  molecular  crose  eections  and 
theoretical  atomic  croes  eectiona. 


(1)  Atomic  cross  sections  calculated  from  Ref.  (I9). 

(2)  Atomic  cross  sections  calculated  by  Hammerling, 

Shine*  astd  Kival  ^19). 
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For  this  calculaition  it  is  neceissary  (and  «ul£icient>  ta  know  th®  comno** 
sitioii  and  density  of  the  ionized  gas.  tho  electron  temperature,  and  the 
average  electron  di£fusion  crocs  cection  oi  all  apecies  which  contribute 
appreciably  to  the  total  electron  collision  probability.  Since  the  results 
indicate  thermodynamic  equilibrium  has  been  achieved  in  those  experi¬ 
ments,  the  temperature,  density  and  composition  of  the  ionized  gas  at 
a  given  shock  velocity  and  initial  air  density,  and  thus  the  concentra¬ 
tions  of  species  which  may  contribute  appreciably  to  the  total  electron 
cross  section  are  known  For  the  molecular  species  and  for  Argon 

Che  electron  diffusion  cross  section  is  available  in  the  literature 
For  atomic  oxygen  and  nitrogen,  the  electron  diffusion  cross  sections 
have  not  been  experimentally  determined,  and  only  somewhat  contra¬ 
dictory  theoretical  estimates  on  the  total  electron  collision  cross  sec¬ 
tions  of  atomic  oxygen  seem  to  be  available  in  the  literature 
Calculations  have  been  made  at  this  laboratory  by  Hammerling,  Kivel 
and  Shine  on  the  electron  collision  cross  sections  of  oxygen  and 

nitrogen  atoms.  The  results  are  consistent  with  the  measurements. 

The  electrical  conductivities  shown  on  Fig.  (21),  particiUarly 
in  the  strong  shock  region,  are  of  an  order  of  magnitude  greater 
than  the  conductivity  of  sea  water.  As  a  resiilt.  the  "thin"  heat  trans¬ 
fer  gage  must  be  insulated  from  the  hot  gas  for  very  strong  shocks  by 
a  thin  layer,  a  fraction  of  a  micron,  of  silicon  monoxide. 

Another  product  of  this  experiment  is  the  indication  that  this 
experiment  may  possibly  be  used  to  determine  the  electron  collision 
cross  sections  of  individual  species  if  the  experiments  are  made  in 
the  component  gases.  It  appears  that  even  in  measurements  in  a 
complex  mixture  such  as  air,  cross  sections  calculated  in  Ref.  (19) 
are  clearly  preferred  over  those  calculated  from  simpler  atomic 
models. 
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MAG  NE  T  OH  YDR  OD  Y  N  AWICS 


Work  in  high  temperatum  gaa  dynamics  and  the  exploitation 
of  the  chock  tube  aa  a  reaearch  tool  natu^-ally  lead  us  to  the  thresh- 
hold  of  magnetohydjrociynamics.  Until  recently  this  field  has  been 
the  private  domain  of  the  aatrophyaicists  and  those  concerned  with 
fusion  reactors >  but  the  interest  of  the  aerodynamicist  in  the  flow  of 
partially  ionized  gaaea  has  injected  new  pointa  of  view. 

The  study  of  the  interaction  between  magnetic  fielda  and 
the  flow  of  electrically  conducting  fluida  covera  a  very  wide  range  of 
physical  phenomena  from  partially  to  completely  ionized  gasea.  At 
the  AVCO  Research  Laboratory  activities  in  both  extremea  of  this 
area  are  in  progress. 

Research  in  the  lower  temperature,  partially  ioniaed  gas 
regirne,  is  aimed  at  exploiting  the  fundamental  advantage  offered  by 
magnetohydrodynamics,  i.  e.,  being  able  to  exert  body  forces  on  a 
fluid  without  any  solid  surface  pushit^on  the  fluid  boundaries,  Magne¬ 
tohydrodynamics,  therefore,  changes  the  basic  flow  equations  eo  as 
to  open  a  new  field  in  fluid  dynamics.  At  the  present  state  of  devel¬ 
opment  of  this  field,  it  is  still  impossible  to  foresee  even  its  most 
important  engineering  applications  (21), 

Experiments  conducted  at  the  AVCO  Research  Laboratory  by 
Patrick,  Brogan  and  Rosa  (^^^  are  designed  to  check  the  one -dimen¬ 
sional  flow  ecuations,  including  the  effect  of  body  forces  due  to  the 
interactions  betv.een  the  conducting  fluid  and  magnetic  fields.  This 
work  can  be  done  in  ordinary  shock  tubes  with  gas  thermally  ionised 
by  strong  shock  waves.  It  is  a  convenient  simplification  to  close  the 
current  loops  in  the  gas,  eliminating  the  electrode  problem  and 
greatly  simplifying  analysis. 

The  principle  of  these  experiments  is  to  create  an  electric 
field  due  to  the  V  x  B  interaction  (velocity  and  magnetic  field)  in  a 
direction  perpendicular  to  both  the  magnetic  field  and  the  How.  In 
the  schematic  of  the  experimental  arrangmert  shown  in  pig.  (22) 
this  field  closes  a  loop  in  the  annulus.  This  current  loop  in  turn 
produces  a  body  force  on  the  flow  in  a  direction  opposing  the  flow. 
The  strength  of  the  magnetic  field  is  arranged  so  that  the  flow  can- 
net  be  decelerated  smoothly  to  sonic  velocity  by  the  body  forces. 

Aa  a  result  a  staztding  shock  must  form  in  the  field  region*  decel¬ 

erate  the  flow  to  subsonic  velocity  and  cause  a  loss  in  total  pressure. 
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When  the  flow  leaves  the  field  region  it  expands  (like  the  driver  of 
a  shock  tube)  but  due  to  the  total  pressure  loss  in  the  standing  shock, 
it  does  not  recover  the  full  velocity  of  the  flow  before  interacting 
with  the  field.  The  change  in  flow  velocity  and  the  location  of  the 
standing  shock  can  be  csdculated  and  compared  with  the  experimental 
results.  A  mirror  camera  picture  of  the  flow  in  the  annulus  under 
the  action  of  the  magnetic  field  i»  shown  in  Fig  (23)  together  with  a 
schematic  explanation  of  the  process. 

At  the  other  extreme,  in  the  fully  ioniaed,  low  density  gas 
regime  eifoTta  are  being  made  to  extrapolate  our  knowledge  to  con¬ 
ditions  interesting  to  the  astrophysicist  and  neuclear  physicist.  Pro¬ 
ducing  new  understanding  in  this  region  requires  first  a  method  of 
producing  in  the  laboratory  a  controlled  sample  of  this  gas  for  study. 
Experiments  performed  at  this  laboratory  to  date  by  Janes  and 
Patschek  have  achieved  temperatures  up  to  400,  000°K  in  the  gas. 
Having  created  such  a  gas  sample,  the  basic  flow  properties  are 
now  being  studied. 

The  problem  of  producing  this  hot  gas  sample  has  been  ap¬ 
proached  through  use  of  essentially  an  electromagnetically  driven 
shock  tube.  A  changing  magnetic  field  is  used  to  accelerate  an 
ionised  mass  of  gas  to  a  high  velocity.  The  stationary  coid  gas  ahead 
of  the  moving  gas  then  has  a  high  relative  kinetic  energy  which  is 
converted  into  thermal  energy  in  a  shock  wave  by  collision  with  the 
moving  gas  and  the  magnetic  field.  Knowing  the  initial  gas  density 
and  shock  velocity,  as  in  the  conventional  shock  tube,  the  conditions 
behind  the  shock  front  can  be  calculated  from  the  conservation  laws. 


One  device  v'hich  has  been  constructed  for  these  experiments 
is  illustrated  schematically  in  Fig.  (24)  and  shown  in  Fig.  (25)  in 
the  illumination  of  its  own  light.  The  discharge  of  a  high  voltage 
capacitor  through  the  primary  coil  causes  current  in  the  secondary 
split  coil  to  produce  an  electrodeless  ring  discharge  in  the  annular 
glass  tube.  This  discharge  is  expelled  down  the  tube  by  the  magnetic 
field.  Our  capacitive  power  supply  operating  at  66  kv  is  capable  of 
10^^  watts.  In  order  to  obtain  temperatures  of  10^  ®K  it  is  necessary 
to  accelerate  deuterium.  to  a  velocity  of  approximately  23  cmj^sec. 

Our  velocity  measurements  to  date  have  indicated  temperatures  of 
up  to  4  X  10^  OK  (14  cm/  ^sec).  l>2  ia  used  because  it  can  be  raised 
to  high  stagnation  temperatures  witlt  a  minimum  of  energy. 
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Fig.  23  —  Mirror  camera  photograph  of  the  luminosity  in  the  flow 
(as  seen  through  the  slit  shown  in  Fig.  (22)  )  for  experi¬ 
ments  conducted  without  and  with  magnetic  field  applied. 
A  schematic  explanation  of  the  important  points  is  also 
shown  for  the  experiineat  with  field  applied.  Note  the 
Steady,  standing  shock  In  the  field  region. 


47- 


PRIMARY  COIL  COMING  FROM 
POWER  SOURCE  (  m  I0'°  WATTS) 


^  SPLIT  SECONDARY  COIL 


- - - AXIS - 


Fig.  24  •—  Simplified  schematic  sketch  of  one  configuration  of  high 
temperature  gas  accelerator.  (Power  delivery  system 
not  shown). 
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Fig.  25  —  Photograph  of  the  existing  experimental  arrangement 
illuminated  by  light  from  the  annular  discharge  tube. 
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CONCLUDING  REMARKS 


The  work  discussed  is  representative  of  some  of  the  work  of 
the  AVCO  Research  Laboratory  during  the  past  two  years.  It  is 
typical  of  many  other  projects  in  which  the  laboratory  is  presently 
engaged  and  demonstrates  the  broad  field  of  application  for  the  shock 
tube  in  high  temperature  gas  dynamics.  It  has  been  shown  that  the 
laroper  environment  can  be  simulated  to  study  many  aspects  of  hyper • 
sonic  flight  problems.  As  an  example  of  an  aerodynamic  problem,  the 
heat  transfer  through  the  laminar  boundary  layer  at  the  .stagnation 
point  has  been  discussed.  A  somewhat  more  physical  phenomena, 
the  emissivity  of  high  temperature  air,  has  been  investigated  and  its 
role  in  hypersonic  flight  has  been  established.  Jn  addition,  the  elec¬ 
trical  properties  of  high  temperature  air  have  been  studied.  In  the 
analysis  of  all  high  temperature  experiments,  the  chemical  kinetics 
problem  has  been  involved  due  to  the  dependence  of  the  flow  param¬ 
eters  on  the  thermochemical  state  of  the  gas.  Experiments  which 
have  defined  this  state  for  the  shock  tube,  and  give  order  of  magni¬ 
tude  answers  for  other  situations,  have  been  described.  In  addition, 
the  preliminary  results  of  our  venture  into  magnetohydrodynamics, 
over  a  large  temperature  range,  ere  presented.  The  scope  of  the 
work  described  shows  the  versatility  of  the  shock  tube  as  a  research 
tool  and  the  need  for  aerodynamicists  to  broaden  their  consideration 
to  include  new  and  previously  unrelated  fields  in  their  attempt  to 
understand  the  problems  of  hypersonic  flight. 
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